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ABSTRACT 

Calculations are presented detailing the gravitational lens diffraction due to 
the steep brightness gradient of the limb of a stellar source. The lensing case 
studied is the fold caustic crossing. The limb diffraction signal greatly exceeds 
that due to the disk as a whole and should be detectable for white dwarf sources 
in our Galaxy and its satellites with existing telescopes. Detection of this diffrac¬ 
tion signal would provide an additional mathematical constraint, reducing the 
degeneracy among models of the lensing geometry. The diffraction pattern pro¬ 
vides pico-arcsecond resolution of the limb prohle. 

Subject headings: gravitational lensing — diffraction — stars: atmospheres — 
white dwarfs 


1. Introduction 

In gravitational lensing two new images appear when a source crosses into a caustic 
curve. The two new images initially appear together on the critical curve with the same 
path length from source to observer. As the source moves across the caustic the images sep¬ 
arate and their respective path lengths change. In microlensing multiple images of a source 
can not be spatially resolved so their amplitudes are combined at the detector producing 
diffraction and interference. The subject of this paper is the diffraction pattern that appears 
near the fold caustic crossing that occasionally occurs during microlensing events. This pat¬ 
tern is proportional to an Airy function, and the Airy diffraction length dj determines the 
periodicity of the intensity oscillations. For sources whose size is less than this characteristic 
length diffraction effects appear at full intensity. For sources larger than this length diffrac¬ 
tion becomes diminished. In this paper we show that for stellar mass main sequence and 
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white dwarf sources significant diffraction fringes will be present in the observed intensity. 
We show detectability of a diffraction signal is dependant on the limb profile rather than 
source size. We apply the results to two observed caustic crossing events, the Galactic bulge 
event 96-BLG-3 and the Small Magellanic Gloud (SMG) event 98-SMG-l using both the 
observed sources and hypothetical white dwarf sources. 

We begin with the point source approximation, describing the lensing geometry and the 
lens equation. Next we describe the time delay surface near a fold type caustic. Using this 
time delay surface we calculate the diffraction integral. We then integrate over extended 
source profiles. We apply these results to various astrophysical situations showing examples 
where the limb effect is significant and we discuss examples where it is not. We compare two 
types of source models, a gaussian brightness profile and a brightness profile that mimics 
the sun. 

Diffraction and interference effects in gravitational lensing have been considered by 
several authors. The theory of gravitational lensing was first derived in a simple manner 
by Refsdal (1964). Some effects of caustic crossings are discussed in Ohanian (1982), 
Blandford & Narayan (1986), and Goodman et. al. (1987). Wave effects and diffraction 
in gravitational lensing are discussed in Deguchi & Watson (1986a) and Deguchi & Watson 
(1986b). Interference is discussed in Peterson & Falk (1991) and Borra (1997). Diffraction 
effects in the frequency domain are discussed in Gould (1992), and Ulmer & Goodman (1995) 
also includes a discussion on white dwarf effects in the Large Magellanic Gloud. Lensing near 
critical points is discussed in Stanek, Paczyhski, & Goodman (1993). Diffraction near fold 
caustics for gaussian sources is discussed in Jaroszyhski & Paczyhski (1995). Determination 
of source limb properties from microlensing is discussed in Gould (2001). 


2. Lensing Geometry 

The gravitational thin-lens geometry for a point source and a point mass lens is shown 
in Figure 1. Although the light is actually in gravitational interaction with the lens over the 
entire path, in the thin-lens approximation there is a single scattering of the light. A ray of 
light moves from the source at point S to the image (or lens) plane. In the plane of the lens 
the ray interacts gravitationally with the lens mass at impact parameter The ray is bent 
by an angle a and travels to the observer at point O. The source position (in the absence of 
a lensing mass) is at an angle (3 from the line of sight to the lens. The observer sees the ray 
at an angle 9 to the line of sight to the lens. All of the angles in the figure are small and we 
assume a Euclidean metric. The distance of the lens from the observer is Dd, the distance 
of the source from the observer is and the distance of the lens from the source is Dds- 
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General relativity predicts that a light ray passing a spherical body of mass M at a 
minimum distance ^ is deflected by the ’’Einstein angle” 

AGM 2Rs 
« = ^ 

provided the impact parameter ^ is much larger than the corresponding Schwarzschild radius 


In the special case of perfect alignment, /? = 0 so the source, lens, and observer positions 
all occur on a line. By symmetry there will then be a ring of images in the lens plane at an 
angle 6 to the line containing the lens and source. This ’’Einstein ring” will be at an impact 
parameter 

Re = = v/2^ (3) 

from the lens using the dehnition D = {DdDds)/Dg. Re is called the ’’Einstein radius”. 

For a point source the ray geometry obeys the lens equation 

f^Dg = ‘^Dds 

Another path, not shown in Figure 1, passes on the other side of the lensing mass 
is described by a second solution to the lens equation. 

In microlensing surveys, lensing light curves are recorded and then an attempt is made 
to determine the mass, position, and relative transverse velocity of the lens. In order to fully 
model a microlensing event one must know all but one of the parameters in equation (4). In 
many cases the lens is not directly observed so the lens position is not known. This leads 
to degenerate solutions of the lens equation and incomplete information on lensing events. 
The lens equation is discussed in detail in Schneider, Fibers, & Falco (1992) and Betters, 
Levine, & Wambsganss (2001). 


( 4 ) 

This path 


2GM 


( 2 ) 


3. Fold Caustics 

If the lens is not a point-like object but consists of several point like objects (i.e. a 
star with planets) the symmetry is usually broken and there can not be a perfect ’’Einstein 
ring”. Instead there can be multiple paths from the source, through the image plane, to the 
observer, each producing an image. If the lens system and source are sufficiently misaligned 
there will only be the image of Figure 1. But as the source and lens move close to alignment 



4 


new images can appear. For a fixed observer and lens, the set of points in the source plane 
where the number of images seen by the observer changes forms a closed curve called a 
caustic. The simplest and most common type of caustic is the fold caustic. This is the 
type of caustic crossing we consider in this paper. The number of images always changes by 
two when a source crosses a fold caustic and the newly created images are generally highly 
magnihed. The newly created images are distinct from the images that existed before the 
source crossed the caustic. 


The set of points in the image plane where the corresponding images are created or de¬ 
stroyed is called a critical curve. Whenever two new images are created they appear together 
on the critical curve. The images initially have the same path length from source to observer 
and therefore interfere constructively. But, as the source moves across the caustic, the im¬ 
ages separate. If the images are unresolved, as in microlensing, the combined image intensity 
oscillates as a diffraction pattern. See Bozza (2000) for a description of the trajectories of 
the images in caustic crossing events. 

It is standard to choose a coordinate system where the fold caustic lies along the y* 
coordinate, the Xg coordinate is orthogonal to the caustic, and x and y describe the image 
coordinates. Near a fold caustic the non-cosmological time delay surface can be approximated 
locally by (Jaroszyhski & Paczyhski 1995) 


r = 


1 


1 3 2 1 2 1 2 ^ 

—x +j/ -x.x-y.y + -x, + ^y,'l 


The quantity 


dc — O' Re 


is the Einstein radius of the lens modihed by the factor 


(5) 

( 6 ) 


a = {K‘^)/{m) 


(7) 


Here K is called the flux factor of the source and (m) is the averaged mass of the lensing 
objects in solar units. The value of a is ~ 1 and for simplicity we adopt that value. More 
precise calculations of alpha requires detailed knowledge of the lensing objects. See Kayser 
& Witt (1989), Witt (1990), and Witt, Kayser, & Refsdal (1993) for details. Equation (5) 
describes the time delay of an image point (x, y) for a given source position {xs,ys)- Choice 
of the origin of the image plane is arbitrary for the purposes of this paper since diffraction 
calculations involve integration over the entire image plane. 


In the geometrical optics limit, the magnihcation of a point source just inside a caustic 
varies as 

^ = ^min + (dc/Xg)^ ( 8 ) 
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where Xg is the distance of the source from the caustic. On the outside of the caustic there is 
only the original image and A = Amin- For a point source on the caustic the magnihcation is 
formally inhnite. However, as seen below, using an extended source and applying diffraction 
theory, the magnihcation remains hnite. 


4. Diffraction 


In scalar diffraction theory the magnihcation of a set of images can be written as the 
modulus squared of the amplitude 


^{xs,yg) = 


^+oo 


^ + 00 




dx / dy expiiojT) 


(9) 


Here the Fresnel length is 

( 10 ) 

The constant in front of the integral normalizes the magnihcation to one in the absence of a 
lens. For a point source near a fold caustic we use the time delay surface given in equation (5) 
in the dihraction integral in equation (9). This approximation only includes the contribution 
of the images that were created in the caustic crossing and so is only valid near the caustic 
crossing event itself. We also adopt a value of = 0 for the distance of the source from the 
caustic in the direction parallel to the caustic. If the size of the source is small compared to 
the radius of curvature of the caustic curve this is a valid assumption. 


Evaluation of the integral leads to 






where Ai is the Airy function and 


d 


dp 




( 11 ) 


( 12 ) 


is the characteristic scale for the argument change, or the Airy diffraction length. Equation 
(11) is a factor of two larger than the result obtained by Jaroszyhski & Paczyhski (1995) 
because we have corrected a minor error in the limits of the integral. This increases the 
magnihcation by four. 
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5. Extended Source Effects 

The physical optics calculation for an extended source requires integration over both the 
source and image planes. Assuming incoherent emission, the contribution to the integrated 
flux of an extended source from the new images near a caustic crossing is given by 

/ +00 ^+00 

dxs dysl {xs,ys)\'^ {xs,ys)f (13) 

■00 J —00 

The integral is over the source plane, and each point in the source plane contributes based 
on an integral in the image plane from equation (9). /(x^, ys) is the surface brightness prohle 
of the source. To proceed further we must choose a model of the surface brightness of the 
source. 

Because the limb brightness gradient will be the dominant factor controlling the amount 
of diffraction, we need a source model with a realistic limb profile. There is little information 
on limb gradients for any star except the sun and a few nearby giant stars. A reasonable 
model radial prohle of the sun can be seen in Figure 2. This prohle of the sun is that described 
in Rhie & Bennett (2002) and Carrol & Ostlie (1980). The intensity is normalized to unity 
at the center. Moving out from the center the intensity falls oh as a cosine. There is an 
abrupt jump to 40% of the central intensity at the stellar radius. 

Assuming spherical symmetry the surface brightness of any point on the surface of the 
sun is then described by 

I{6) r 0.4 + 0.6 cos(6*) r < Tq 
Jo \ 0 r > r© 

The angle 6 is related to the distance r from the center of the star as seen on the projected 
solar disk by 

9(r) = sin~^— (15) 

re 

This model is a good approximation to the observed solar prohle. We can use the 
discontinuous approximation because for all cases in this paper the width of the region over 
which the limb intensity increases is much less than the dihraction length df (Lites 1983), 
(Kuhn at. ah 1998). 

When considering white dwarf sources we scale this solar prohle to the smaller radius of 
the white dwarf. This assumes that the strong gravitational held at the white dwarf surface 
maintains a steep density gradient in the white dwarf atmosphere. 


There are two scales for dihraction ehects when a solar-prohle source crosses a caustic. 
The hrst is the disk dihraction ehect of the entire source as it moves inside the caustic region. 
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Just inside the caustic, when the source radius is much smaller than df, all components of 
the source are localized within one fringe of the Airy function. As the source moves further 
into the caustic the entire source contributes nearly identical values of the Airy function and 
there is full amplitude diffraction. When a source is larger than df different regions of the 
source will be at different phases within the Airy function. This leads to an averaging, so 
for large sources diffraction effects become diminished. 

The second effect is a limb effect. As the source first begins to cross the caustic only 
its limb is magnihed. This is shown in Figure 3. Most of the source has not crossed the 
caustic and those points are not yet additionally imaged. So when a source of any size 
hrst crosses a caustic its limb will always show diffraction effects. However, even for the 
limb of the source, the averaging effect will begin to take place as the limb moves into the 
caustic. The diffraction oscillations will begin to average out when the limb has passed many 
characteristic lengths df into the caustic. 

In order to extract the small limb diffraction signal from the large unlensed contribution 
a differential observation is appropriate. One method is to rapidly measure differences of 
intensity as the source moves across the caustic. By sampling the event repeatedly such that 
for every time step the source has moved a distance Xg < df, and then subtracting successive 
intensity measurements, diffraction due to a portion of the source radius less than df can be 
isolated. We have calculated this differential magnihcation for several cases described in the 
next section. 


6. Results 

To begin, we compare our calculations to previous results. Fold caustic diffraction 
effects for gaussian sources in the quasar lensing system Q2237+0305 have been studied by 
Jaroszyhski & Paczyhski (1995). In this event, a source at redshift Zg = 1.695 is assumed 
to be lensed by solar mass objects at Zd = 0.0394. The system has a = 0.8. At wavelength 
A = 500nm, df = 4.0xl0®cm. Using these values and a gaussian radial prohle (scaled to the 
appropriate source size), we evaluate the integral in equation (13) and obtain the results in 
the left panel of Figure 4. The magnihcations we calculate are larger than those previously 
calculated due to the difference in the evaluation of the diffraction integral (equation [11]) 
and proper normalization of the gaussian (Jaroszyhski, M., and Paczyhski, B. 2002, private 
communication). For sources with a characteristic width greater than df diffraction effects 
are not visible. 

To study the effects of other source prohles on diffraction we replace the gaussian prohle 



with a model solar profile. Although the source in the case under consideration is a QSO, 
which does not have a solar prohle, we use this prohle to show that steep brightness prohle 
gradients enhance diffraction. The right panel in Figure 4 is for solar prohle sources of 
various radii. For sources larger than dj diffraction fringes have a greater amplitude for a 
solar prohle source than a gaussian prohle, particularly at large Xs- For sources with radii 
of lOdf, for either prohle, the dihraction ehects are not visible. For these larger sources we 
need to look at the diherential magnihcation to see fringes. 

Limb dihraction is minimal for gaussian prohle sources but is much stronger for a solar 
prohle source. Figure 5 shows the magnihcation and diherential magnihcation of a gaussian 
source with a characteristic width a = lOdf as it passes into a fold caustic. There are no 
apparent dihraction ehects in Figure 5. Figure 6 shows the magnihcation and diherential 
magnihcation of a solar prohle source with a much larger radius (lOOOd/) as it crosses a 
fold caustic. Figure 6 shows dihraction ehects in the diherential magnihcation plot. For 
sources 1000 times the Airy dihraction length, even though the whole-disk dihraction ehect 
is almost completely averaged out, a substantial limb ehect can be seen. The steep gradient 
of the scaled solar prohle produces the dihraction signal seen in Figure 6 . Because a gaussian 
prohle has much smaller gradients, or alternatively because the width of the limb region is 
larger than df, the limb ehects are severely diminished. 

There are several ongoing microlensing surveys using sources in our Galaxy and its 
satellites. These surveys occasionally report caustic crossing events. The MACHO/GMAN 
collaboration (Alcock et. ah 2000), the PLANET collaboration (Albrow et. ah 2001), the 
OGLE experiment (Udalski et. ah 1992, 1997), and the MPS collaboration (Rhie et. al. 
1998) are all currently using sources in either the Galactic bulge or Magellanic clouds. We 
now adopt lensing parameters applicable to these surveys. 

The event 96-BLG-3 was the third bulge caustic crossing event detected by the MAGHO 
project in 1996. A I.IRq GO star at 8.5kpc was lensed by a I. 6 IM 0 binary lens at 7.1kpc. 
The velocity of the source perpendicular to the caustic was 103kms~^. This event is described 
in Alcock et. al. (2000) and Lennon et. al. (1996). The event 98-SMG-l was the hrst 
SMG caustic crossing event detected by the MAGHO project in 1998. The time of the 
second caustic crossing (passing out of the multiple image region) was predicted far enough 
ahead of time that all of the groups listed above were able to observe this event. A I.IR 0 
star at 62.5/cpc was lensed by a 0.36Mq binary lens at 58kpc. The velocity of the source 
perpendicular to the caustic was 76kms~^. This event is described in Alcock et. al. (2000) 
and Rhie et. al. (1999). We use these two observed events as our canonical examples. 

First we take the 96-BLG-3 case. Evaluating the diffraction integral we get the results 
in Figure 7. Since d/ is the typical length for diffraction effects this is the resolution obtained 
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by measuring these effects. For this system df = 9.6 x lO^cm, or 7.5 pico-arcsecond on the 
sky. For the quoted source velocity of 103kms~^ the initial period of oscillation is 1/lOs. 
The disk magnihcation reaches a maximum of ~ 78. Although this is slightly less than 
the observed value of ~ 120 (Alcock et. ah 2000) the source passed near a cusp and 
therefore the peak magnihcation would be slightly higher than we predict when assuming a 
fold crossing. As seen in the differential magnihcation plot, the amplitude of the dihraction 
fringes is ~ 5 x 10“^ of the unlensed source intensity. 

Next we take the 98-SMC-l case. Evaluating the dihraction integral we get the results 
in Figure 8. This system has d/ = 2.3 x lO^cm and translates to 2.49 pico-arcsecond on the 
sky. For the quoted source velocity of 76kms~^ the initial period of oscillation is l/3s. The 
disk magnihcation reaches a maximum of ~ 78. This is in good agreement with the best ht 
to the data of ~ 70 (Rhie et. al. 1999). As seen in the diherential magnihcation plot, the 
amplitude of the dihraction fringes is ~ 1 x 10“^ of the unlensed source intensity. 

The signals produced in the 96-BLG-3 and 98-SMC-l cases are too small to be detected 
with available technology. Since compact objects produce strong dihraction ehects we now 
use white dwarfs as the sources in our models. We replace the main sequence sources in 
our two previous cases with a IMq DA white dwarf similar to Sirius B, with a radius of 
0.0087?©. Using this white dwarf as the source object in the 96-BLG-3 event we evaluate 
the dihraction integral and get the results in Figure 9. The data points were separated by 
source displacement steps of size d//4. The disk magnihcation reaches a maximum of ~ 800. 
In the diherential plot, the amplitude of the dihraction fringes is ~ 8% of the unlensed 
source intensity. Substitution of a white dwarf source does not change the initial period of 
oscillation or df. 

Finally we take the case where the same white dwarf is the source object in the 98- 
SMG-1 lensing case. Evaluating the dihraction integral we get the results in Figure 10. The 
disk magnihcation reaches a maximum of ~ 1100. As seen in the diherential magnihcation 
plot, the amplitude of the dihraction fringes is ~ 40% of the unlensed source intensity. 


7. Discussion 

For the cases we have considered the whole disk dihraction signals are undetectable. For 
main sequence stars in the Galactic bulge and SMG the limb dihraction intensity is ~ 10“^ 
the unlensed source intensity, also undetectable with current technology. However, for white 
dwarf sources in the Galactic bulge and SMG limb dihraction is 10 to 40% of the unlensed 
source intensity. It is interesting that dihraction is more readily detected in small sources 



than in large ones. This is opposite the usual bias in astronomical observations. 

So far there have not been any reported white dwarf caustic crossing events. The re¬ 
ported events have much brighter source stars. However, if white dwarf events are discovered 
in the Galactic bulge or Magellanic clouds it should be possible to detect a diffraction signal. 

We have considered only a single wavelength in this analysis. However, the effects we 
describe may be easier to detect if broad optical bands can be used. All the diffraction 
patterns presented here scale with wavelength so if a broad bandwidth is to be used it must 
be subdivided into spectra. The resolving power required is approximately the number of 
diffraction fringes from the caustic to the limb of the source. Only low resolution (A/AA ~10) 
is needed. However, these spectra must be recorded in rapid succession since the fringe 
pattern changes in a fraction of a second. 

Caustic crossing events are usually discovered after the source crosses into the caustic, 
which does not allow time to set up rapid spectrographic observations. Precise measurement 
of the event, along with modelling, allows one to predict when the source will pass back out 
of the caustic. The diffraction pattern produced as a source leaves a caustic is reversed in 
time compared to the pattern produced as it enters. The inbound crossing can be used to 
alert observers to record rapidly sampled spectra during the outbound crossing. 

Detection of limb diffraction would provide a measurement of df resulting in an addi¬ 
tional constraint which helps remove the degeneracy from the lens model. Also, accurate 
measurement would allow fine detail of the limb brightness profile to be studied. Diffrac¬ 
tion patterns could provide information on distant collapsed objects with higher angular 
resolution than is available through any other technique. 
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Dd Dds 


Fig. 1.— Point source, thin-gravitational-lens geometry. The source is at point S, 
the lens at point L, and the observer at point O. The ray is bent by an angle a at an impact 
parameter ^ from the lens. The observed image appears at an angle 6 to the line of sight to 
the lens. The source position is at an angle (3 to the line of sight to the lens. Another path, 
not shown in the figure, passes on the other side of the lens in the same plane. 
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Fig. 2.— Model brightness profile of a solar mass main seqnence star. The intensity 
falls as a cosine until the limb is reached at which point it falls from 40% of the central 
intensity to 0. 





Fig. 3.— A circular source crossing a fold caustic. The source coordinates Xg and Ps 
are shown. The caustic is the solid horizontal line and the source is moving up, into the 
interior of the caustic. The source positions are shown at three different times. The source 
moves a distance Xg between positions. The areas A and B that move into the caustic at 
each step are shown. At position 2 area A produces the diffraction signal. At position 3 area 
B produces the dominant signal. 
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Fig. 4.— Magnification of various gaussian profile (left) and solar-profile (right) 
sources. The magnification is plotted versus the distance the source has passed into the fold 
caustic. The geometry is that for Q2237+0305. The plots are for sources with characteristic 
widths[radii] of O.l(solid), 0.3(dotted), l(dashed), 3(dot dashed), and lOdf (dot dot dashed). 
For the solar prohle, the source distance Xg is measured from the limb as in Figure 3. For 
the gaussian prohle, the source distance Xc is measured from the center. In both cases the 
distances are scaled by the Airy diffraction length df. For a gaussian prohle, for source 
widths comparable to df, dihraction is strongly attenuated, particularly for large Xc- For a 
solar prohle source of size comparable to df fringes are still visible at large Xg. In both cases, 
for sources with radii less than the characteristic size df, dihraction ehects are visible. In 
both cases, for sources of size lOdf, no dihraction is visible in these plots of magnihcation. 




Fig. 5.— Magnification and differential magnification for a gaussian source. The 

source width is lOdf. The lensing geometry is that of Q2237+0305. The distance from the 
source center to the caustic Xc is scaled by the Airy dihraction length df. The right panel 
shows diherences of amplitude taken at intervals df/A. No dihraction fringes are visible. 
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Fig. 6.— Magnification and differential magnification of a solar profile source. 

The source has a radius of lOOOd/. The geometry is that of Q2237+0305. Xg is the distance 
the limb of the source has passed into the caustic. In the magnihcation plot this distance 
is scaled by the radius, in the differential plot it is scaled by the Airy diffraction length df. 
The corresponding scale expansion is rsource/df=lQQQ. The right panel shows differences of 
amplitude taken at intervals d//4. Even for a source of radius 1000 times the diffraction 
length, limb diffraction effects can be seen for this solar prohle source. 




Fig. 7.— Magnification and differential magnification for the event 96-BLG-3. 

Xs is the distance the limb of the source has passed into the caustic. In the magnification 
plot this distance is scaled by the radius, in the differential plot it is scaled by the Airy 
diffraction length df. The corresponding scale expansion is rsource/df=7.99 x 10"^. The plots 
are sampled at d//4. The peak-to-peak amplitude of the diffraction signal in the differential 
plot is ~ 5 X 10“^ of the unlensed source intensity. 
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Fig. 8.— Magnification and differential magnification for the event 98-SMC-l. 

Xs is the distance the limb of the source has passed into the caustic. In the magnihcation 
plot this distance is scaled by the radius, in the differential plot it is scaled by the Airy 
diffraction length dj. The corresponding scale expansion is rsource/df=3.29 x 10"^. The plots 
are sampled at df/4. The peak-to-peak amplitude of the diffraction signal in the differential 
plot is ~ 1 X 10“"^ of the unlensed source intensity. 




Fig. 9.— Magnification and differential magnification for a hypothetical white 
dwarf source in the 96-BLG-3 event. Xg is the distance the limb of the source has 
passed into the caustic. In the magnihcation plot this distance is scaled by the radius, in 
the differential plot it is scaled by the Airy diffraction length df. The corresponding scale 
expansion is rsource/df=581. The right panel shows differences of amplitude taken at intervals 
df/4. The peak-to-peak amplitude of the diffraction signal in the differential plot is ~ 8% 
of the unlensed source intensity. 
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Fig. 10.— Magnification and differential magnification for a hypothetical white 
dwarf source in the 98-SMC-l event. Xg is the distance the limb of the source has 
passed into the caustic. In the magnihcation plot this distance is scaled by the radius, in 
the differential plot it is scaled by the Airy diffraction length df. The corresponding scale 
expansion is rsource/df=24Q. The plots are sampled at d//4. The peak-to-peak amplitude of 
the diffraction signal in the differential plot is ~ 40% of the unlensed source intensity. 





